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Abstract 

The simulation methods for the channehng phenomena in GeV/TeV energy range 
>^ I in ideal or distorted crystal lattices are discussed. Monte Carlo predictions for feed- 

out and feed-in rates, dislocation dechanneling, and deflection efficiencies of bent 
Oh. crystals are compared to the experimental data. The role of multiple interactions 

(— I I with crystal in circular accelerators (" multipass channeling" ) for the efficiency boost 

in the crystal-aided extraction experiments is analysed. Possible future applications 
of the crystal channeling technique are considered. 



1 Introduction 

The experiments on charge particle beam steering by means of channeling in bent 
crystals have greatly progressed in recent years, spanning over two decades in energy[||, 
^, ^. The detailed theory of beam steering is essentially based on Monte Carlo sim- 
ulations, as tracking of a particle through a bent crystal lattice requires not only a 
calculation of a particle dynamics in this nonlinear field, but also a generation of 
random events of scattering on the crystal electrons and nuclei. 

To track particles through the curved crystal lattices in simulation we apply the 
approach with a continuous potential introduced by Lindhard. The continuum-model 
description is disturbed by scattering which could cause the trapped particle to come 
to a free state (feed out), and an initially free particle to be trapped (feed in). By 
every step the probabilities of scattering events on electrons and nuclei are computed 
depending on their local densities which are functions of coordinates. This ensures 
correct orientational dependence of all the processes in crystal material. 
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2 Numerical methods 



As the channeled particles move far from the nuclei, the transitions between the 
random and channeled states are mainly due to electronic scattering. It is these 
processes which have been studied experimentally in detail. Hence, an adequate 
theory needs an accurate consideration of electronic scattering Q|. 
The energy transfer T in a collision has a distribution function: 

<fN _ Dpe{x) 1 



dTdz 2/52 T2 ■ 

Here D=0.3071 MeV civ?g~^ Zf Z / A; [3=v/c; the local density pe{x) of electrons is 
normalized on the amorphous one. Deviations from (||) at r~/ or T^iT^ax are of no 
concern, due to the nature of dechanneling as explained below. The (round) angular 
kick from collision is 

X 1/2 

2meT{l-T/T„,ax)] Ip (2) 



Tmax — lTae(? 0^^"^ is the maximal transfer to a single electron. The frequent small 
kicks produce a diffusion-like angular scattering, with mean square 

^Ims = ^(Yl ~ Ti/Tmax)) (3) 

^ i 

Relativistic Diffusion Factor 

The diffusion approximation for electronic scattering was used by all the authors (e.g. 

P). One assumes that scattering on electrons is diffusion- like, i.e. 6s <C 9c in any 
collision. Then the angle of particle is changed by small steps, with the rms value 
from Eq.(P). Note that 9rms depends on the total transferred energy, not on the detail 
of (H). In the MeV range this approximation for heavy ions works well, because even 
the maximal possible angular kick per collision 9'^'^^ = y/nieTmax/p ~ lArUe/M {M 
for the particle mass) is < 9c- 

At ~100 GeV (the range of modern applications of bent crystals), 9c is greatly 
reduced to ~10 /irad, so rare catastrophic collisions with 9s>6c may happen. The 
problem for the diffusion approach is that the integration up to Tmax in the diffusion 
factor (^) is no longer justified. The energy transferred with catastrophic collisions 
{9s>9c) is of no importance for dechanneling, and therefore it should not be included 
into Eq. (|3|). Although 9^ms depends on T^ax via ln{Tmax/I), the removal of the 
energy transfers of catastrophic collisions from Eq. (^) reduces the diffusion coefficient 
by a factor of ~2. Hence, the diffusion approach may overestimate the rate of feed- 
out and feed-in by about the same factor. Here is an example. For a 100 GeV proton 
Tmax =10 GeV. However the transfer Tc causing the angular kick (projection) equal 
to 9c, is quite moderate: 

^ p202 2M7^ 
rric rric 

Tc= 4 MeV for a 100 GeV proton in Si. From the energy transferred via single 
collisions, about one half is due to the scatterings with 9s>9c- For collisions with 
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T>Tc the T value is no longer important. E.g., scatterings with T=10 MeV and T=l 
GeV are equally essential as they knock the particle out of channeling mode at once. 

The above algorithm, i.e. consideration of electronic scattering as a series of sin- 
gle scattering events instead of the traditional diffusion approach, is the pronounced 
feature of the computer code CATCH|Q]. Other approaches, such as analytical cal- 
culations by Gartner K. et al.Q, can nicely describe the effects dominated by the 
potential of crystalline lattice, notably the lattice distorted by defects, but the mat- 
ter of our study - the transitions from channeled to random states and backward in 
perfect bent crystals as induced by electronic scattering, - can be approached most 
correctly by Monte Carlo methods. 

The predictions of CATCH were testified by the recent experiments, where feeding- 
out, feeding-in, high-efficiency bending, energy loss spectra, dislocation dechanneling, 
and crystal extraction from accelerators have been studied [Q. 

3 Crystal extraction: multi-pass channeling 

In the extraction mode, circulating particles may pass through the crystal many 
times. Another major point: the first incidence of particles may occur very close to 
the crystal edge. The scattering of unchanneled beam in crystal, and the accelerator 
optics become important. 

3.1 The SPS Experiments 

Before the CERN SPS studies of crystal extraction Q, theoretical comparisons @ 
with extraction experiments |l^ were restricted by analytical estimates only, which 
gave the right order of magnitude. The computer simulations considered idealized 
models only and predicted the extraction efficiencies always in the order of 90-99% 
while real experiments handled much smaller efficiencies, in the order of 0.01%. 

The consider ed-below theoretical work has been the first comparison between 
the realistic calculation from the "first principles" (computer simulation) and the 
experiment. The simulation was performed [0] with parameters matching those of 
the SPS experiment. Over 10^ protons have been tracked both in the crystal and 
in the accelerator for many subsequent passes and turns until they were lost either 
at the aperture or in interaction with crystal nuclei. In the simulation two options 
were considered. The first one assumed near-surface irregularities (a 'septum width') 
of a few fim. With impact parameter below 1 /im, it excluded the possibility of 
channeling in the first pass through the crystal. The second option assumed perfect 
crystal surface. Table |l| shows the expected extraction efficiencies for both options 
from the first simulation run and the measured lower limit of extraction efficiency as 
presented at the 19-th meeting on "SPS Crystal Extraction" Q held at CERN. 

Though the efficiency comparison, theory to measurements, was not possible at 
that time, from the analysis one could see that the perfect-surface simulation pre- 
dicted narrow peaks for the angular scans (30 ^urad fwhm) and extracted-beam pro- 
files, which were not observed. The imperfect-surface option was consistent with the 
experiment: wide (~200 //rad fwhm) angular scan and sophisticated profiles of the 
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Table 1: SPS crystal extraction efficiencies from the early runs, Monte Carlo and experiment 



Option 



Monte Carlo Experiment 



Poor surface 



Ideal surface 



15% 



40% 



lower limit 

of 2-3% 
only known 



extracted beam (dependent on the crystal alignment). The efficiency was then mea- 
sured in the SPS with that first tested crystal to be lOzbl.7%. Detailed simulations 
have shown that efficiency should be a function of beam vertical coordinate at the 
crystal, and be from 12 to 18% at peak, with imperfect-surface option. 

The simulation studies for a new crystal with "U-shaped" geometry, performed 
prior to the measurements, predicted just slight increase in efficiency and much nar- 
rower angular scan in the option of edge imperfection. Fig. |l] shows one of measured 
scans in rather good agreement with prediction. Further simulation took more realis- 



tic details into model [13|, and studied the dependence of crystal efficiency on septum 
width. Table |2|. These simulations have been repeated with the energies of 14 and 270 
GeV, where new measurements have been done at the SPS. The results are shown in 
Table H vs measured data. 



Table 2: Peak efficiency F vs septum width t. The statistical error is 0.6%. 



t (/im) 



20 50 100 200 



F 



13.9 12.4 12.9 10.9 8.2 



3.2 The Tevatron Experiment 

The Tevatron extraction experiment has provided another check of theory at a sub- 
stantially higher energy of 900 GeV. A detailed report of predictions for this ex- 
periment from the Monte Carlo simulations was published^ two years before the 
measurements were taken |2[. 

We have investigated three options: a crystal with ideal surface, one with a septum 
width (amorphous layer) of t=l /im, and one with t=50 /xm. The crystal bending 
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Figure 1: The angular scan of SPS extraction with a U-shaped crystal. Prediction ((8>) 
and measurement (:*r). 

shape and other details were as used later in the experiment. Figure ^ shows that 
there is little difference between the three options; the peak efficiency is about 35- 
40%, and the angular scan fwhm is 50-55 /urad. This insensitivity to the crystal 
surface quality is due to the set-up different from that used in other experiments; as 
a result, the starting divergence of incident protons at the crystal was not small and 
hence less sensitive to edge scattering. 

The measured peak efficiency was about 30%. This value, together with the 
measured angular scan, is superimposed in Figure ^ on the theoretical expectation, 
showing a rather good agreement. 

3.3 Crystal optimisation 

The length of the Si crystals used in the SPS and Tevatron experiments was about 
optimal to bend protons with a single pass. The efficiency of the multi-pass extraction 
is defined by the processes of channeling, scattering, and nuclear interaction in the 
crystal, which depend essentially on the crystal length L. The length optimisation 
was a subject of a simulation, with results shown in Fig.^. 
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Figure 2: Angular scan of the efficiency at Tevatron. Ideal crystal (•); imperfect crystal: 
septum width t=l fim {-k), t=50 /im (*). Also shown is the measured peak efficiency and 
angular scan. 



3.4 Protvino crystal extraction experiment 

This experiment, started at the end of 1997, aims to test crystals as short as possible 
in order to gain in the extraction efficiency from an increased number of proton 
encounters with crystal, as promised by simulations. In the first run of 1998 this 
experiment established world record of efficiency of crystal extraction, over 40% §. 
The silicon crystal was just 5 mm along the beam direction, and only 3 mm long 
central part was bent to give the 70 GeV protons a deflection of 1.5 mrad. Figs. 
show the extraction efficiency and the angular scan, Monte Carlo vs measurement. 



4 Future Applications 

To understand further possible applications of the technique it is useful to recall the 
analytical theory of multipass crystal extraction [0]. Few simple assumptions were 
taken as an input to the theory: any particle always crosses the full crystal length; 
pass 1 is like through an amorphous matter but any further pass is like through a 
crystalline matter; that there are no aperture restrictions; particles interact only with 
the crystal not a holder. 

The overall multipass channeling efficiency is then derived to be 

where 

^{L/Ln) = exp{-kL/LN) ~ (ttLn/L)^/^ - 1.5 (6) 

may be called a "multiplicity factor" as it just tells how much the single-pass efficiency 
is amplified in multipasses. A fraction (l-G) of the channeled particles is to be lost 
along the bent crystal due to scattering processes and centripetal effects. Then the 
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Figure 3: Extraction efficiency vs crystal length as simulated with imperfect surface. 
Tevatron (top) and SPS (bottom). Also shown is the measured efficiency at 4 cm. 



multipass extraction efficiency is 
Fe = FcxG 



X T.{L/Ln) X G 



(7) 



, 2 J (Jidp 

The theory check against the CERN SPS data |15] shows good agreement (Table H). 



Table 3: Extraction efficiencies (%) from the SPS experiment, theory [17], and detailed 
simulations. 

pv{Ge\) SPS Theory Monte Carlo 



14 0.55±0.30 0.30 0.35±0.07 
120 15.1±1.2 13.5 13.9±0.6 
270 18.6±2.7 17.6 17.8±0.6 



From (^) we see that multiplicity factor can be huge if L is very small ox big. 

MeV extraction. One opportunity (small L) is inspired by the recent successful 
experiment [^] on bending 3-MeV proton beam by means of graded composition 
Sii-xGcx/Si strained layers. This invention allows to cover the whole spectrum of 
accelerator energies (from MeV to multi-TeV) by bent crystal channeling technique. 
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Figure 4: Efficiency of Protvino crystal extraction as a function of the fraction of beam 
store incident on the crystal; the measurements (crosses) and simulations (open circles). 



F(%) 




-0.5 



Angle (mrad) 



Figure 5: Angular scan of Protvino crystal extraction: prediction (•) and measurement 
(o). 



8 



Now one can consider extraction from accelerators starting with MeV energies, by 
crystals as short as from 1 /im. Eqs.(||-P) predict that channeling efficiency over 99% 
can be achieved in sub-GeV (and up to several GeV) range, thus opening a new 
world for bent crystals applications. With traditional bent crystals, it was common 
to think that highest efficiencies are achievable at highest (TeV) energies as multiple 
scattering angles vanish with energy faster than channeling angle does. It is very 
interesting now to find that channeling efficiency can be even more boosted at lower 
energies due to huge multiplicity factors. One can build a very efficient system to 
extract beams from accelerators with crystals. 



Muon extraction. The other opportunity (big Lj\f) for huge multiplicity factor 
is muons which have formally = oo. Theory then says that efficiency of muon 
channeling should be 100%. Actually the multiplicity factor for muons is limited by 
muon lifetime (mostly) and muon outscattering of the aperture. With a muon mean 
lifetime of 1000 turns[|l^, the number of encounters with crystal may be big. 

At muon machine, the backgrounds "have the potential of killing the concept of 
the muon collider" 1 17 ] . As muons cannot be absorbed, it was proposed to extract 2- 
TeV muon beam halo with electrostatic septum as a primary element. Positive muons 
can be easiely steered away by bent channeling crystals. Short analysis says that we 
could steer negative muons also, e.g. by the same bent planes as used for positive 
particles, Si(110)||l8[]. In same crystal Si(llO), dechanneling length is shortened 
by factor of ~100 for negative particles. However, at 2 TeV is huge (~1 meter) 
for positives and modest (~1 cm) for negatives. The required deflection angle is only 
64 //rad[16] and can be ensured by a Si crystal ~1 mm long — quite shorter than L^. 
Again, multiplicity factor greatly favors muons of both sign. One can build a very 
efficient system to handle halos at muon colliders. 
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